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Abstract 

Two component superlattice film of 2N layers is considered and dimensionally quan- 
tised spectrum of phonon-s is found. 

The problem of localization of phonon-s in the superlattice with random thicknesses 
of the layers is investigated. The Landauer resistance of the transport of phonon-s is cal- 
culated exactly. For short range disorder the numerical analyses shows, that at frequency 
u=0, there is delocalized state and the correlation length index v is equal to 2. 
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1 Introduction 



The interest to superlattice structure in condensed matter physics is based on the possi- 
bility of manipulation of the physical properties of devices by changing of characteristic 
lattice parameters. 

The growth techniques can be used to prepare specimens consisting of alternating 
layers of thickness d\ of constituent A and thickness d<i of constituent B. Samples can 
be prepared so that d\ and o?2 have any value from two or three atomic spacings up to 
hundred atomic spacings. The entities A and B can be materials with different acoustic, 
electronic, magnetic properties and one can consider semiconductor, metal, insulator and 
superconductor constituents, attempting to change values of expectable physical variables 
in a desirable regime. 

The technical advance in fabrication of superlattices motivated an intensive study of 
various physical properties of these systems, especially electronic and vibration spectra, 
optical and magnetic properties, e.t.c. The calculations P, || shows that the spectrum 
of quasi-periodic systems are intermediate between periodic and random ones. 

Along with electronic properties study of elastic waves in a bulk superlattices has 
been a subject of interest in the last decades |18|, [19|, |2l], |22|, ^3|, |24|, |25|, |27|, [28|, p9fl . 
A one dimensional theory of acoustic vibrations in layered material was given long ago 
by Rytov M. Elastic waves have also been investigated in semi-infinite superlattices 



m be m, 33, 33 



The consideration of the superlattice films instead of the massive ones provides ad- 
ditional opportunities for controlling the elastic and electronic parameters of the super- 
lattices. A sufficiently complete experimental knowledge about the oscillator spectra 
of binary laminated semiconductor systems is available nowadays p4 |. Numerous and 
generally mutually compatible results about the frequencies of long-wave phonon-s in 
InAs — GaSb, Ge — GaAs are presented in literature . 

Together with strongly periodic superlattices the effects of localization and tunnelling 
of the electrons was studied in short rang disorder superlattices [|, ^j, [7| 
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The problem of localization of electrons in random potential and hopping parameters 
in low dimensional spaces are in continuous interest of physicists after Andersons remark- 
able article ||. Originally Mott and Twose conjectured that all states are localized in 
ID systems for any degree of disorder. It was argued [T(J that in case of full randomness 
all states are localized in dimensions equal and less of two. However, recent investigations 



shows [JTT], [12], [13] that delocalized states can appear in case of correlated disorder. 

The aim of present article is twofold. First we found the spectrum of transversal 
phonon-s in two component superlattice film with boundaries with arbitrary finite number 
of slices. The results, obtained here, are also applicable for longitudinal waves, which 
moves in the transversal to the layers direction. Second, we consider random distribution 
of the thicknesses of the superlattice compounds and calculate the Landauer resistance 



14|1 of the acoustic phonon-s exactly. The continuous model is used, in which the layers 



are considered as a macroscopic elastic bodies. 



2 Notations and equation of motion in the bulk 



We follow here the notations and derivations of the book [pqj . 

Let the vector of elastic displacement (deformation) of the material at the space 
point x is u(x). For small variations (u << 1) we can derive the strain tensor (tensor of 
deformation) as 

The response of the free energy F to variation of deformation tensor defines the stress 
tensor as 

dF = —SdT + Uikdiiik, (2) 

According to Hooke's law, for small deformations mk around equilibrium, we can restrict 
ourself to consider only quadratic in terms in series expansion of F. Define 

1 



dF 



Kuiia ik + 2(JL [ - -UiAk 



dui k , 



(3) 
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which defines the following stress tensor 

a ik = KuuCTik + 2/i (iiik - -SikUi^j , (4) 

where K is the modulus of compression and fi is the modulus of rigidity. K > and 
[i > 0. One can express stress aik in terms of Young's modulus E = and Poisson's 
coefficient a = \ as follows 



E 



a 



Oik = z— t — u ik + - —u u 5 ik . (5) 

1 + a \ 1 — la J 

Then, the equation of motion for the displacement «j is simply 

da^ . . 

P u i= ^ — 6 

OXk 

where F$ = is the force, which acts on unite volume around x with density of matter 

p(x). After substituting (|J) into (|6|) we can define the equation for Ui(x) in a bulk. 

Ej E Bui 

P u i= zAui H — : — rgrad-— — . 7 

H 2(1 + a) 2(l + a)(l-2a) & dx l y ' 

For the transversal waves 

divu = 0, (8) 

and from the fl7|) we can easily get 

d 2 Ui 



dt 2 



c 2 t A Ui = 0, (9) 



with 



For the longitudinal waves 



E 



< = w^ry (10) 



Au = grad(divu) (11) 
the equation (0) reduces to the following equation 

^ - c 2 e An, = (12) 

where q is the longitudinal velocity 

E(l-a 



p(l + a)(l-2a)' 



(13) 



3 Boundary conditions in the superlattice film, Trans- 
fer Matrix and the spectrum 



Let's consider superlattice, elementary cells of which consists of two layers of various 
materials A and B with the thickness d 2 and modulus of rigidity fii, fi 2 (fig-1)- The 
number of pairs in the film is N. 
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Fig.l. 



We consider transversal elastic waves (diva = 0) propagating inside superlattice in 
arbitrary direction. It can be shown, that all results are reproducible for longitudinal 
waves if they propagates in transversal to layers direction. 

Lets choose coordinate system such, that waves are propagating in the (x, y) plane (x 
represents transversal to layers direction), with the wave vectors (k\, q±, 0) and (k 2 , q 2 , 0) 
in the A and B materials correspondingly. Without loose of generality one can take 

U X =Uy = 0. 

The solutions of the wave equations (^) with frequency u, which full-fills transversality 
condition divu = 0, is the superposition of forward - and a backward-travelling waves 

«2n-i = (c 2 n-ie iklX + c^e'^) e^-^, 

u 2n = (c 2n e ik2X + c 2n e- ik2X ) e< qy - wt \ n = 1 • • • N, (14) 

where for hi, i = 1,2 we have 

k* + q 2 = 4 (15) 

c it 
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In formulas (0) 2n(correspondingly2n — 1) numerates the layers of B(or A) type and 

Cit{c2t) are the velocities of sound in that materials. 

We should now impose the boundary conditions on the displacements u 2n and u 2n -x- 
Let's consider now free boundaries of the film, which means the use of Neumann 

boundary conditions 

d x ui = d x u 2 N = 0. (16) 

On the boundary of the A and B layers one should use the continuity condition for 
the displacements 

u 2n = u 2n -i, (17) 

as well for the forces 

Fi = a ik ds k . (18) 



In the equation (|18D ds k = n k ds-is the normal vector to the boundary and equal to small 
area in modulo. Hence we have 

°ik n k = °ik n k- (19) 
On the boundaries of the film the forces are equal to zero 

Vikn k = (20) 



By use of expression (|) for a ik from the boundary conditions (|l6l)-|T7|) and fll9| - |2"0"D one 
can easily obtain following set of equations for the displacements u 



ul N = < = 0, 
H 2 uf = inul n -\ 



u 2n = u 2n-\ n =l,2,---N. (21) 

This equations transforms into the following equations for the coefficients of the forward 
and backward travelling waves 



C 2 N 



e ik 2 (d 1 +d 2 )n _ ^ -ifc 2 (di+d 2 )n _ g 



fi 2 k 2 C2ne ik2[idl+d2)in - 1]+dl] - fx 2 k 2 C2ne- ik2[idin+d2{n - 1)] = 
= ^ l k l c 2n -ie ikl[din+d2{n ~ l)] - nihc 2n -ie~ ikl[din+d2( - n ~ l)] 

ik 2 [d 1 n+d 2 (n-l)] , -= -ik 2 [d 1 n+d 2 (n-l)] _ 
L 2n c ~r <^2n c — 

n „ik 1 [d 1 n+d 2 (n— 1)] _i_ ^ _— ik\[d\n+d2{n— 1)] 

L 2n— l c ~r <^2n— l c 



ci — ci = 



We will solve this set of linear equations by use of transfer matrix method [|18| 
Let's define now 



In 



('2, 



\ C 2 r 



Then the half of the set of equations (p3|) can be reformulated as follows 



Mn^2n = B 2n -1^2n-l- 1 



with 



.4 



2n 



gik 2 (ndi+(n—l)d 2 ) ^,—ik 2 (ndi+(n—l)d 2 ) 

^ik 2 {nd-i+(n-l)d 2 ) e —ik 2 (ndi+(n-l)d 2 ) 



and 



B 2n -i 



Mi fei gtfci (ndi +(n— l)d 2 ) g— ik-j (nrfj +(n— 1)^2) 

Ai2fe ' M2&2 



D jfci (ndi + (n— 1)1^2) 



-jfci (din+(n— 1)^2) 



Equation (24) can be rewritten as 



1p 2n = A 2 l B 2n -! 1p2n~l- 



Similarly, the other half of the equations (^) looks as 



■02n-l = ^2n-l B 2n -2 1p2n-2, 



with 



.4 



/ „ik 1 ((n-l)di+(n-l)d 2 ) 



2n-l 



^—iki((n-l)di+(n-l)d 2 ) 



gifci((n— l)di + (n— 1)^2) g— ifci((n— l)di+(n— l)d 2 ) 
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and 



2n-2 



M2fc 2 gifea +(n-l)rf 2 ] 
Ml fc l 

_ifc 2 [(n— l)di+(n— l)d 2 ] 



M2fc 2 e — ifc 2 [(n-l)d] +(n-l)rf 2 ] 
-ifc 2 [(n-l)di+(n-l)<Z 2 ] 



Recursion equations (P7|) and (ESI) allows us to connect ?/> 2n with ^1 in a following 



^ 2n = A 2 l B 2n _ x A 2 l_ x B 2n _ 2 ■ ■ ■ A 2 1 B a Vi- 



Let's now define the Transfer Matrices as 
?i = B 2n _i A 2 ^_ 1 



/ cos Mi sin Mi 

/J2&2 M2«2 



and 



^2 — -B 2n An,, 1 



1 i sin k\di cos fcirfi 



^cosM 2 sin k 2 d 2 

z sin /c 2 c? 2 cos fc 2 c? 2 



Then the equation (31) becomes 



ij>z n = B£T n A 1 il> 1 = Uil> 1 



where 



cos fcic?i cos k 2 d 2 — i cos fcidi sin k 2 d 2 + ^ 

— ia^i sin fcidi sin A; 2 ci 2 , sin fcidi cos fc 2 <i 2 



T = T X T 2 



i sin fcidi cos k 2 d 2 + 



^2/02 



sin fcidi sin k 2 d 2 + 



\ +^^" cos ^1^1 s i n ^2^2 + cos fclG?l cos fc 2 c? 2 
The first equation of fl22|) for n = N can be written as 



where 



^ 2Ar = (g*2(iVdi+./Vd 2 ) _ e -ift 2 (iVdi+7Vd 2 )^ 



For the another boundary of the superlattice film, where n = 1 we have 



which means that 



1 



(39) 



From the equations Q3"3) ) and Q3D| ) we can obtain following equation for the spectrum of 
transversal phonon-s in the superlattice film of 2N layers 



Tr 



CT 



A' 



0. 



(40) 



where 



1\Q 



(41) 



To proceed further we need to calculate the N — 1 degree of the Transfer Matrix, 
which can be achieved simply by diagonalizing T. Obviously 



A JV , 



0, A 



N 



w, 



(42) 



where A and A are eigenvalues of T, and W is the diagonalizing matrix. One easily can 
find the eigenvalues or the Transfer Matrix T as 



A = e 



cos kidi cos k 2 d 2 — - \ — -r- + ^—7- I s i n k\d\ sin k 2 d 2 ) ± 



2 V fx 2 h Hik\ 



± i 



\ 



1 — ( cos k\d\ cos ^2^2 — - f — \ 1 + 2 ) sin /^g^ sin /c 2 ^2 ) , (43) 



2 \/X2&2 jUi^i 



where 



cos O = cos k\d\ cos k 2 d 2 — ( ^ 1 1 + ^ 2 2 ) sin &iG?i sin k 2 d 2 . 

\ji 2 k 2 HikxJ 



(44) 



Further, a simple calculations shows, that equation (f40"D reduces to 



/mA^ = 0, 



(45) 



which means that 



1 • ■ • N. 



(46) 
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Finally we obtain following equation for the spectrum of transversal phonon-s 

cos k\di cos k 2 d 2 — | ^ , 1 + ^ 2 , 2 ) sin k\d\ sin k 2 d 2 = ± cos n^- (47) 
V/X2K2 HikiJ ' ' N 



where 



k-=K- q 2 (48) 



We see that this equation is coinciding with the equation for the spectrum of phonon- 
s in the bulk [I], 0> but the momentums in perpendicular to the layers direction are 
quantised due to dimensional restriction of the film. 

4 The Landauer resistance of phonon-s in the super- 
lattice with random distribution of thicknesses of 
the layers 

The problem of elastic waves in superlattice is essentially one dimensional. One di- 
mensional problems are especially attractive because of their possible exact integrability. 



In the article by Erdos and Herndon |L5| the problem of the transport of particles in the 
one dimensional space for a wide class of disorders was considered in the Transfer Matrix 
approach and general results were obtained. It was proved that that Transfer Matrix of 
the one dimensional problem belongs to SL(2, R) group and randomness can be exactly 
taken into account for such quantities as Landauer resistance [14]. 

Some exact results for Kronig- Penney model in case of non-diagonal disorder by other 



methods was obtained in fli"7|. 



It is easy to see from the formulas ([35]) for the Transfer Matrix T, here we also have 
a representative of the SL(2, R) group. One can make a link between transfer Matrices 
of the Kronig-Penney model and phonon-s in the superlattice. 

Following jnj and |I5| let's define resistance as a ratio of reflection to transmission 



coefficients, which, by use of formula (0), is equal to 



|2 



P = = U l2 U* 12 = U l 2 (U + )l (49) 



T 
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where JJ\ is the 1,2 matrix element of the evolution matrix U 

U = B^(hT 2 ) N Ai. (50) 

We are going to consider random distribution of thicknesses of the layers and take 
the average of Landauer resistance. For further convenience we will normalise Ti(T 2 ) 
Transfer Matrices on order to have a unit determinant. It will not change the equation 
(p0|) because the normalisation factors for T\ and T 2 cancels each other. Hence we will 
consider 

' (^t) 1/2 cosA; 2 (x 2i -x 2i _ 1 ) i(^yi/2)smk 2 (x 2l -x 2i ^) 
v *(^|) 1/2 sin k 2 {x 2l - x 2i _0 (g£)VS» cos k 2 (x 2t - x 2i _0 

for the even slices. The similar expression for odd slices T 2i _i can be found simply by 
permuting variables k and fi for 1 and 2. 

Now let us analyse the direct product of the evolution matrices, the (U®U + ) 2 '1 matrix 
element of which defines Landauer resistance. For this purpose we should calculate first 
the simplest constituent block of that expression, namely direct product Tj ® of Ui-s. 



T, 



2i 



(51) 



In the article |T5j it was demonstrated, that this direct product can be represented as 
1 © (3x3) =4 x 4 matrix. It happened because of the fact, that Tj matrices are a spinor 
representations of the SL(2, R), hence, the direct product of two 1/2-representations can 
be expanded as a sum of scalar and vector representations. In the language of the group 
elements T G SL(2, R) this expansion looks like 

cwcrr 1 )? = - ^)2Af>% (52) 



where 



Af = -TrpfpTr 1 ^) (53) 



is the spin-one part of the direct product. But for Landauer resistance we need to 
calculate T ® T + . It is easy to see from the formula ( |51~D that 



g{T- Vx = T+, (54) 
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therefore, by multiplying the the expression fl52|) in the left and right by o\ we will have 

1, 1 



(55) 



Now the calculation of the direct product U ® U + is straightforward. The product 
Ili=L Ti of Tj-s transforms into product of A^-s. Finally we will obtain 



N 



2(^0y(ii A *-i A *r(^i: 



i=l 



(56) 



Substituting this expression, together with the expressions for B 2 ^ and A\ (from (j30|) 
and (|29| ) correspondingly) into the (p£9|), after some simple algebra for Landauer resistance 
p we will have 



P 



1 M\ki 



2 fi 2 k 2 
+ i(A 



-1 + (A 



A? 



^1 1 + ^2^2 

2 p 2 k 2 pih 



) + 



AT\12 1 /t^lkl P2k 2 + 



(57) 



2 /i 2 &2 P\ki 

where (A^) 11 (correspondingly (A^) 12 ) is the 11(12) matrix element of the matrix A = 
A1A2, which is a product of A-s of the I and II slices. 

The average over any type of random distributions of the parameters of the model 
can be calculated now exactly. We consider random distribution of thicknesses of the 
slices, keeping boundaries fixed xo = 0, x 2 n = L. We see from the formula ( 51 ) that Tj 
depends only on the thickness of the slice x-i — x^\. The only restriction we have is the 
condition, that 

2N 

E A ^ = L - (58) 
i=i 

Therefore, the average of the A^, with the probability distribution g(y) ; (/ °° g(y)dy = 
1), defined in the following way 

N ,„ 2N N 

/ 



N 2N N 

(II A 2 i-iA 2i > = / dy 1 . . .dy 2N g{yi) . . .#(t/ 2 aO<KX)% ~ L ) II A 2i-i(?/2i-i) A 2i (?/ 2i ) 
i=i Jo j=i i=i 



^e-^((A 1 (p))(A 2 (p)» 



A' 



(59) 
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where 

(AiM) = / dye^g(y)A 1>2 (y). 
Jo 

The average Landauer resistance is equal now 



(60) 



(P) 



2 \i 2 k 2 



_1 _|_ ij'Pi^i _|_ ^ 2 ^ 2 ' 
2 fi 2 k 2 fiih' 



dpe^ L ((A 1 (p))(A 2 (p))) 



JV 



ii 



+ 



1,/il/Ci /i 2 /c 2 . 



dpe** «Ai(p)>(A 2 (p)» 



iV 
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(61) 



2 /i2^2 Pl^l' 

It is obvious, that in a case of homogeneous media (two components of the superlattice are 
coinciding) we restore the expression for the Landauer resistance of electrons, obtained 
in 
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For large sample size (N » 1), as it was argued in [14 , |16[| , the resistance should 
behave as e N ^, where £ is the correlation length. Excitations are localized or not depends 
on the behaviour of £. If at some frequencies correlation length becomes infinite, we have 
delocalized state and the expression (|61]) shows, that the answer depends on average 
value of A^. For further analyse let's consider simplest case of the distribution, namely 
when there is equal probability for slices to have a thickness up to di (i=l,2). 



g{y) 



di ' 
0. 



< y < di 
otherwise 



(62) 



We have taken Al x Gai_ x As and GaAs as a components of the superlattice with the 



parameters [34| 



li x = (3.25 - 0.09x)10 u dyn/cm 2 , \i 2 = 3.2510 u dyn/cm 2 , 
pi = (5.3176 — 1.6x)g/cm 3 , p 2 = 5.3176(?/cm 3 , 

d x = 30 ■ (5.6532 + 0.0078x)A°, d 2 = 10 ■ 5.6532A 



(63) 



and consider waves, propagating in the perpendicular to the layers direction (<f = 0). 

For large enough iV the asymptotics of p, and therefore the correlation length £, 
defined by the closest to unity eigenvalues of (Ai)(A 2 ). If it is A, then 



~ l/lnA(cj). 
12 



(64) 




Fig.2 

This result is easy to understand, uo = means that we have constant displacement u, 
which simply is the sift of the all sample. The correlation index v from £ ~ , defined 
as a slop of the plot of lnln\{uj) versus Inuj is presented in Fig.2 and it appeared that 
u—2. All other states are localized. 
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